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Abstract 

We examine the published hght curves (LCs) of 117 X-ray binary pulsars, focusing on the dependence of 
their light curves on the observed energy bands. It is found that the energy dependence of the LCs appears 
only when the X-ray luminosity is larger than ~ 5 x lO'^^erg s~^. Assuming that the behavior of light 
curve is related to the radiative accretion column on the neutron star surface, this energy threshold can be 
considered as the observational proof of the accretion column formation proposed by Basko and Sunyaev. 
Once we can grasp the existence of radiative column, we can also obtain several useful informations on the 
neutron star properties. As an instance, we perform the statistical analysis of the orientation angle of the 
magnetic axis, and we find that the inclination angle of magnetic axis should be small in order to explain 
the observed statistics. 
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1. Introduction 

After the discovery of the pulsating X-ray sources 
(Giacconi et al. 1971; Tananbaum et al. 1972), a huge 
number of observational and theoretical studies of X-ray 
binary pulsars (XBPs) have been published. The basic 
concept of X-ray pulsars has been understood in a quite 
early stage of these studies; that is, the periodical X-ray 
pulsation is considered to be originated from accreting 
magnetized neutron stars (theoretically, Pringle & Rees 
1972; Lamb et al. 1973; observationally, White et al. 
1983; Nagase 1989; Bildsten et al. 1997). 

In the standard model of X-ray pulsars, the X-ray 
source is a rotating neutron star (NS) with a binary com- 
panion which feeds accretion matter via Roche lobe over- 
flow or stellar wind. If the NS has a strong magnetic 
field (^ 10^^ G), the accreting matter eventually cannot 
accrete onto NS surface across the magnetic field lines. 
Rather accreting matter is trapped to the magnetic field 
and constrained to follow the field lines, and near the stel- 
lar surface, it forms a funnel-like flow onto the magnetic 
poles (e.g. Pringle & Rees 1972, Basko & Sunyaev 1976). 
Consequently, the accreting matter is concentrated near 
the polar caps and emits high-energy radiation as the X- 
ray photons. Since generally the magnetic axis and the 
rotational axis are misaligned, these X-ray photons from 
the poles are observed as periodic pulsations (Shapiro & 
Teukolsky 1983; White et al. 1983; Nagase 1989; Bildsten 
et al. 1997). 

The basic model of X-ray pulsars has already been 
obtained. Although the detailed investigations of X-ray 
sources in these systems are still difficult since the obser- 
vational method of XBPs is limited in space experiments. 



the sample number of XBPs increase significantly in these 
days, thanks to recent development of space experiments. 
Hence, even though what we can know from a single ob- 
servation is less, we can distill some hints of X-ray source 
properties by statistical analysis of observed X-ray pul- 
sars (Bulik et al. 2003; Raguzova & Popov 2005). For 
instance, Bulik et al. (2003) use the catalogue of XBP 
and propose the new open question about the distribu- 
tion of peaks of X-ray LCs concerning with the formation 
scenario of the NS binary systems. They have counted 
the number of single and double peak X-ray LCs in pub- 
lished catalogues, and found that the observed number of 
single peak pulsars is relatively high. If the inclination 
of magnetic axis from the rotational axis is random, the 
majority of observed LC should have double peaks. Then, 
they argued that, in order to explain the large number of 
single peak LCs, the inclination angle of the magnetic axis 
with respect to the rotational axis must be small. 

In the investigation of the NS physics, LCs of XBPs give 
us useful informations of the systems. The unique strong 
point of LC observation is that it does not depend on the 
distance to the source and carries us original information 
of central NS. Generally, the light curves of those XBPs 
are various. Some of them show simple single/double peak 
pulses persistently, while a lot of systems shows compli- 
cated light curves. These complexity is caused by absorp- 
tion by circumstellar matter, misalignment of the mag- 
netic axis and rotation axis, Doppler effect, and so on. 
Also the shape of the LCs often depend on the observed 
energy band (White et al. 1983; Nagase 1989). Several 
mechanisms have been proposed in order to explain the 
energy dependency of LCs. For example, if the matter 
surrounding the NS has absorption rate which depends of 
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the energy band, the observed X-ray LC will be different 
in each energy band (White et al. f983; Marsden et al. 
1998). However, the real situation is still under discussion. 

In this paper, we examine whether wc can say some- 
thing new about NS physics from recently provided XBP 
catalogues. We mainly focus on the LC of XBP, which 
brings us the distance-independent information about the 
magnetized NS. Firstly, we check the relationship of LC 
shapes with other observational properties of XBP. We 
find there is clear correlation between the LC complexity 
and the X-ray luminosity of XBP. This correlation can be 
explained by classical radiation column model given by 
Basko, Sunyaev (1976). Then wc make statistical anal- 
ysis of the LC of XBP. Statistics of the LC distribution 
can be considered to reflect the axis inclination conditions 
of central NS. Our analysis suggests that the inclination 
angle between the magnetic axis and the rotational axis 
of the central magnetized NS should be small. Also it in- 
dicates that the opening angles of the X-ray beams from 
the radiative column should be narrow. Since these NS 
properties arc difficult to obtain from direct observations, 
such a statistical investigation becomes important. 

In the next section, we introduce the X-ray binary pul- 
sar samples briefly. In Sec. 3, we remind the basic model 
of radiative zone of NSs. In Sec. 4, the data and analysis 
of X-ray pulsar samples is shown. In Sec. 5, we analyze 
the XBP sample statistically and discuss the obtained na- 
ture of central magnetized NS. The last section is devoted 
to the summary. 

2. Analysed systems of XBP 

Thanks to several space observations by ASCA, Bcppo 
SAX, CGRO, Chandra, Ginga, INTEGRAL, ROSAT, 
XMM-Ncwton and so on, now wc have a large number 
of XBP samples and candidates. Firstly, we should men- 
tion that we are interested in the X-ray pulsars whose 
LCs are published. These XBPs are listed in several cata- 
logues; here we refer to the catalogue of Be/X binaries 
and candidates given by Raguzova, Popov (2005), the 
catalogue of high-mass X-ray binaries given by Liu, van 
Paradijs, van den Heuvel (2000), the review and tables of 
discrete X-ray sources in Small Magellanic Cloud (SMC) 
given by Yokogawa et al. (2003), the table of Be/X-ray 
binaries in SMC given by Haberl, Pietsch (2004), the 2nd 
IBIS/ISGRI soft gamma-ray survey catalogue given by 
Bird et al. (2005), and the list of observed light curves of 
X-ray pulsars given by Bulik et al. (2003). From these 
previous works, we chose 117 X-ray pulsars whose spin 
periods are known. 70 of them are Galactic sources, 41 
are in SMC and 6 are LMC sources. These chosen sys- 
tems are listed in Tables. 1-3. Even comparing with the 
generic and comprehensive previous analysis of X-ray pul- 
sars given by White et al. (1983), Nagase (1989), Bildsten 
et al. (1997), Bulik et al. (2003) and so on, our sample 
number is the largest. 

Wc prepare the following three tables; systems in 
Galaxy, Small Magellanic Cloud (SMC), and Large 
Magellanic Cloud (LMC). In these tables, the X-ray bi- 



nary systems arc sorted by the spin period. In the third 
column of the tables, we show the type of light curves, 
which is categorized as following. Among all the sam- 
ple XBP, for 90 systems, the LCs have been taken and 
published in at least 2 energy bands. Then we catego- 
rize those pulsars into 2 groups: Regular LC group and 
Irregular LC group. The diagnosis is the following. When 
we compare two light curves obtained in different energy 
bands, if (i) the number of peaks arc different, and/or (ii) 
the phase of the peak is shifted (roughly more than 20 
%), then we categorize this pulsar as Irregular group ^. 
Otherwise, it is put into Regular group. Some pulsars 
show marginal behaviors of LC, hence firstly we catego- 
rize them into Regular group but indicate them as "Rm" 
in the Tables. For the systems whose light curve is pub- 
lished only for one energy band, we can not make diagnosis 
and indicate " — " in this column. 

The fourth and fifth columns are devoted to the spin 
period of the NS and the orbital period of the system, 
respectively. The spin period Pg is shown in unit of sec- 
ond, while the orbital period Pq is in unit of day. In 
the sixth column, the X-ray luminosity Lx (erg s~^) is 
shown. Note that many systems of XBP are observed as 
transient X-ray sources. For almost all the systems, Lx 
coincides to the maximum luminosity of the system, and 
their LCs are taken in this bright phase. The surface mag- 
netic field strength of the NS Bp, which can be estimated 
from cyclotron line, is shown in seventh column. In the fi- 
nal column, we show the additional information. Also the 
references of light curves are indicated here. We have to 
keep in mind that those data has been taken by different 
experiments, and hence the total luminosity of these XBPs 
are estimated by different energy bands. This difference 
may make uncertainty of the total luminosity shown in 
the tables. 

3. Radiative column 

The luminosity and mass accretion rate is strongly 
connected in NS accreting systems. In Basko, Sunyaev 
(1976), they derived the "accretion rate" as 

L,.M^. (1) 

where M is mass-accretion rate. And they proposed a 
critical accretion rate which is described as 
„ 2cloGM 
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^ For dim sources and sources without enough statistics, small 
peaks of LCs could be missed. Hence there may bo a tendency 
that dim XBPs are easy to be counted as Regular group. On 
the other hand, the diagnosis of the peak shift has no bias with 
regard to the luminosity. 
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In the first equation, k denotes the opacity. Here, Cs de- 
notes the isotropic scattering cross-section, and ctt is the 
Thomson cross-section, respectively. indicates the cir- 
cumference of the accretion cylinder (see, Fig. 1 in Basko, 
Sunyaev 1976). They predict that if Lt is smaller than L*, 
then the accretion channel is wholly optically thin. Hence, 
the X-ray radiation is mainly emitted from the small cross 
section of the accretion column on the surface and the 
main beam is emitted into the direction along to magnetic 
axis. On the other hand, if Lt is larger than L*, the op- 
tical depth of the accretion column above the NS surface 
will be large. This optically thick region is called radiative 
column. In this radiative column, the radiative pressure 
is no longer negligible and the accreting fluid decelerates 
significantly. Passing this radiative column, the fluid con- 
verts its kinetic energy into radiative energy. Since the 
emitted photons cannot go straight in this dense region, 
the photons are scattered and diffusively go out from side 
face of the column. Here, note that, since the geometry of 
the accretion flow is not spherical, the radiation pressure 
becomes important even when the luminosity is less than 
the Eddington limit for spherical accretions. 

The emission manner depends on the existence of ac- 
cretion column. If the opaque radiative column appears, 
the X-ray radiation will be emitted from the side-face of 
the column, while the radiation is emitted in the polar 
direction from the cross-section of the column when such 
a column does not exist. X-ray emission in the direc- 
tion of magnetic pole is called the pencil-beam, while side- 
face radiation is called the fan-beam. Considering these 
two manners of X-ray radiation, we can expect that the 
critical accretion rate L* can be understood also as the 
threshold of exchange between these two types of radi- 
ation. Actually, there are examples which shows the de- 
pendence of the LC on the X-ray luminosity. For instance, 
the time dependent LC and Lx of EXO2030-h375 (G28) 
has been taken and the results show that the phase of the 
primary peak is completely changed with the threshold of 
Lx « 103'^[erg s'^] (Parmar et al. 1989a,b) This shows 
the configuration of the accreting column changes at this 
threshold luminosity and subsequently the beam type is 
changed (Nagasc 1989). 

4. Data analysis of XBPs 

First of all, we check whether the observational quan- 
tities of XBP show any difference between Irregular and 
Regular groups. In Fig. 1, the luminosity distributions 
of Irregular group, Regular group, and the whole sam- 
ple of XBP are shown by thick solid, thick dashed, and 
thin dotted lines, respectively. From this figure, we can 
see that Regular and Irregular groups have different dis- 
tributions (with more than 99% probability according to 
the Kolmogorov-Smirnov test). The Regular group shows 
almost symmetric distribution with an axis of symme- 
try around Lx ~ 2 x 10^^ [erg s~^], and the upper and 
lower tails spread broadly. On the other hand, Irregular 
group distributes in a narrow region above Lx > 5 x 
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Fig. 1. Number counts of Regular and Irregular pul- 
sar luminosities. The luminosity distributions of Irregular 
group, Regular group, and total sample including pul- 
sars whose LC is unpublished, are shown by thick solid, 
thick dashed, and thin dotted lines, respectively. Prom 
this figure, we can see that the distribution of Regular 
and Irregular groups are different. Two exceptionally 
dim Irregular pulsars are G30 {XTEJ1906-I-09) and S28 
(RXJ0059.3-7223), which are considered to be wind-fed type. 

Irregular LC pulsars, and they are G30 (XTEJ1906+09) 
and S28 (RXJ0059.3-7223). Both of them have giant com- 
panions and are fed by wind accretion. 

Here we note that, especially for Galactic sources, it is 
difficult to estimate the distance of the source precisely. 
Because of this uncertainty, the luminosities of Galactic 
XBPs may have errors as large as factor of 10 (Haberl 
& Sasaki 2000). On the other hand, the distances (and, 
hence, luminosities) of SMC/LMC XBPs can be estimate 
much precisely. In order to avoid the uncertainties of the 
luminosities of Galactic XBPs, we can check the number 
counts of Regular and Irregular XBP luminosities only 
for SMC/LMC sources, and confirm that the resultant 
distribution is almost the same with Fig. 1. Also we would 
like to comment that the uncertainties of the luminosity 
could broaden the distributions not only of Regular group 
but also of Irregular group in Fig. 1. Consequently the 
inconsistency of the distributions will be still essential, 
even with the uncertainty. 

In Fig. 2, the relationship between the spin period Pspin 
and the luminosity is shown for Regular (open-circle) and 
Irregular (fiUed-box) LC pulsars. For the purpose of the 
comparison, also No-LC group pulsars are shown by small 
cross in the same figure. Almost every pulsar distributes 
on upper-right region, while a few exists in lower-left re- 
gion. These exceptional points are corresponding to the 
pulsars which arc in the quiescent phase or in the propeller 
regime (Bildsten ct al. 1998). According to the analysis 
given by Ikhsanov (2003) and Raguzova, Popov (2005), 
the accretion mode exchanges from the accretion to the 
propeller regime on the basis of a certain criterion. That 
is, if the condition 
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Fig. 2. Spin period distribution regarding to luminosity is 
shown for Irregular (star) and Regular (cross) LC pulsars. 
For comparison, No-LC group pulsars are also shown (plus). 
Almost every pulsar distributes on upper-right region, while 
a few exists in lower-left region (they are in the propeller 
regime, or in the quiescent phase). The critical line of the 
propeller transition, for fi = 10^" is also shown. It indi- 
cates that in the propeller regime, the LC must be regular. 
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Fig. 3. Spin period of Irregular and Regular LC pulsars 
are plotted as the function of the magnetic field strength 
on the stellar surface. No significant difference can be 
seen between Irregular and Regular groups. The excep- 
tionally strong magnetic field pulsar with ~ 10'^*[G] is G14 
(1E1048. 1-5937) and it is perhaps an AXP. We may see 
weak dependence of spin period on the magnetic field. 



will be prevented. Here, usq and L^q are the normal- 
ized magnetic moment and luminosity, respectively. At 
a glance, the Irregular pulsars distribute almost in the 
same region as Regular pulsars. Notice that, however, 
entire Irregular group locates above, or near, the criti- 
cal line of the propeller transition, shown by Eq. (3). It 
indicates that in the propeller regime, the LC becomes 
regular, and it is consistent with the argument of Basko, 
Sunyaev (1976). Namely, in propeller regime, the accre- 
tion rate must be much smaller, since the accreting matter 
will be blown away. In this regime, and with consequent 
small accretion rate, it seems difficult to form the radiative 
dense column. Hence, in this case, the X-ray radiation is 
mainly emitted as a pencil beam from tiny region on the 
stellar surface. Then, in such regime, we can only observe 
simple pulse emitted in pencil manner. 

Fig. 3 shows the spin period of Regular and Irregular LC 
pulsars plotted regarding to the magnetic field strength on 
the stellar surface. No significant difference can be seen 
between Regular and Irregular groups. Since X-ray pulsar 
systems are not equilibrium systems, they are considered 
to be still in spin-evolutionary stage. It means that the 
spin period needs not any relation with magnetic field. 
There is one point which has exceptionally strong mag- 
netic field (~ 10i''[G]), and this is G14 (1E1048. 1-5937). 
In our sample, G13, G14, G68 and S9 are the candidates of 
anomalous X-ray pulsars (AXPs). If they are, the radia- 
tive mechanism could be different from XBPs, hence we 
should remove them from our sampling. In the analysis, 
we do not have the necessary LCs for G13 and S9, and 
those pulsars are not included in the statistical results. 
G68 is a reliable candidate of AXP, hence we exclude this 
from all analysis. On the other hand, G14 is counted as a 
sample, hence this one point can be an error. Historically, 
G14 has been considered as a candidate of XBP in the 
early stage (Seward et al. 1986), however, it is supposed 
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Fig. 4. Number counts of spin periods of Irregular 
and Regular pulsars. The spin period distributions of 
Irregular group, Regular group, and total sample includ- 
ing pulsars whose LC is unpublished, are shown by thick 
solid, thick dashed, and thin dotted lines, respectively. 

as a candidate of an AXP in the last few years (Tiengo et 
al. 2002). 

In Fig. 4, the spin period distributions of Regular group. 
Irregular group, and total number including pulsars whose 
LC is unpublished, are shown by thick solid, dashed, 
and dotted lines, respectively. The number of observed 
Regular pulsars increases with spin period and truncates 
at around 1,000 sec, while, the number of Irregular group 
shows rather flat distribution. The spin of the magne- 
tized NS accelerates via accretion process and decelerates 
in propeller regime. The Irregular group pulsars, however, 
locate in the parameter region where propeller effect does 
not work (see. Fig. 2). Hence, those Irregular pulsars do 
not suffer from the spin down mechanism. Then, they may 
keep their spin rapid and show rather flat distribution as 
shown in Fig. 4. 
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Fig. 5. Number counts of orbital periods of Irregular and 
Regular pulsars. Lines are the same as in Fig. 4. The orbital 
period number-counts of Irregular group shows quite flat dis- 
tribution, while Regular group has a peak around ~ 100[d]. 

In Fig. 5, tlie orbital period distributions of Regular 
group, Irregular group, and the total number including 
pulsars which LC is unpublished, are shown by solid, 
dashed, and dotted lines, respectively. From this figure, it 
is clear that Irregular group shows quite flat distribution 
of the orbital period, while Regular group has somehow a 
peak around 100 d. Unfortunately, it is difficult to observe 
long orbital periods, say, longer than 1,000 d, since they 
need a long observational time (several years). Hence, 
near the upper end, we may have observational bias. Here, 
we may expect the following scenario: the longer orbital 
period means the larger orbital separation. The larger 
separation may cause smaller accretion rate, if the mass 
of the primaries are almost the same. Hence, from the 
condition given by < L*, the number of Irregular group 
decreases in large Porbit region. Although this expectation 
matches to the tendency of the figure, it is required better 
statistics to confirm it. 

5. Discussions 

As shown in Fig. 1, Irregular group has clear minimum 
X-ray luminosities around Lx « 5 x lO'^^erg s^^, and this 
minimum limit coincides with the criterion of radiative 
column given by Basko, Sunyaev (1976). Namely, it is 
natural to consider that the irregularity of LC correlates 
with the existence of radiative column with finite height. 

Although there are two exceptions in Fig. 1, these two 
XBPs have different type of companions and different ac- 
cretion manners: wind-fed accretions from large compan- 
ions. That is, the central NS of these systems are sur- 
rounded by huge amount of circumstellar matter. Because 
of the large absorption by the surrounding matter, the 
small peak may be screened. Also the dumpiness of the 
surrounding matter may cause the phase-dependent con- 
cealments of LCs (Nagase 1989). These effects will be 
significant in lower energy bands, and they lead to the 
energy-dependent LC shapes. Therefore, it may be pos- 
sible to conclude that the irregularity of these pulsars 



is caused not by the configuration of radiative column, 
but by absorptions due to surrounding matter. If this is 
the case, all the other Irregular pulsars satisfy the Basko- 
Sunyaev criterion, and can be considered to have radiative 
column with finite height. 

Examining strictly, we have to treat the X-ray lumi- 
nosity carefully. The X-ray luminosities taken by differ- 
ent experiments are estimated based on different energy 
bands. For XBP, it is difficult to estimate the total lu- 
minosity accurately, only from fragments of the spectra. 
However, in our sample, almost Irregular XBPs are ob- 
served in relatively broad energy bands, from 0.1 keV to 
several tens keV. It means that the X-ray luminosities 
of Irregular XBPs are estimated based on large part of 
the total spectra, with relatively high accuracy. Hence, 
our conclusion that the Irregular group has large X-ray 
luminosity is unchanged in our sample, even under the 
consideration of the uncertainty of the luminosity due to 
the lack of the spectra. 

Reminded that, if a NS does not have a radiative col- 
umn. X-rays will be emitted only as pencil beams in our 
basic conjecture. Otherwise, X-ray LC is affected from 
the contribution of fan-beams. Next we check that the ra- 
diative column really affects the shape of LC. And then, 
we examine whether we can obtain further informations 
of magnetized NS from the statistical analysis of LC. 

5.1. Light curve models 

The X-rays from a NS are radiated in two different man- 
ners, that is, pencil beam and fan beam. The intensities of 
these two beams have the angular dependence; the pen- 
cil beam gives the maximum intensity when the line of 
sight coincides with the magnetic axis, while the fan beam 
gives the maximum intensity when we observe the NS from 
the magnetic equator. Hence, as the simplified models of 
beams, we assume that the intensity (number of photons) 
can be described as 

,x f cos^a (pencil beam) ,.^ 
^ ^ 1 sm a (tan beam), ^ ' 

where a denotes the angle between the line of sight and the 
magnetic axis, fi and v denote certain factors generally 
depending on the spin phase ^. This cos - sin dependence 
comes from the radiating area seen by the observer (Nagel 
1981). We assume simple dipole magnetic poles, and do 
not include the effect of screening by the stellar surface, 
hence the both beams are observed at — 7r/2<Q;<7r/2. 

In general. X-ray photons with only higher energies can 
be emitted as the pencil beam, since along the magnetic 
axis the mean free path becomes so short that weak ra- 
diation can not escape in this sense, when the accretion 
column exists. On the other hand, from the side face 
of the accretion column, photons scattered and diffused 
many times are emitted as the fan beam. These side-way 
photons can have lower energies in average. When the 
height of the column is significant, the fan beam will be 
dominant, however, with a moderate height of column, 
pencil beam could not be neglected. Hence, the observed 
beams of X-ray from the accretion column will be the mix- 
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ture of the fan beam and the pencil beam. We model such 
a mixed beam as following;. 

F{a,ip)^Acos''a + Bsm''a, (5) 

where A and B are mixing factors which are defined by 
A + B = 1 and A G [0, 1]. Of course, ^ = corresponds 
to the fan-beam, while B = results in the pencil-beam. 
In Fig. 6, we show some examples of the light curves gen- 
erated in this manner, with different Om and 6o- Here, 
^ni is the inclination angle between the magnetic axis and 
the rotational axis, while 9o is that between the rotational 
axis and the line of sight. 

If the radiative region has finite side-face, it is natural 
to consider that both beams arc emitted in some amounts 
(mixed-beam) . When we observe the mixture of them, the 
mixed beam will show more complicated pulse shapes. If 
we admit that the ratio A/B depends on the energy of 
photons, the shape of the LC will be changed for each 
energy band, and the energy dependence of X-ray LCs 
can be explained naturally. For instance, if we see Fig. 6 
, and compare pencil-, fan- and mixed-beams, it is eas- 
ily seen that some of them show the phase inversion and 
peak splitting. Such a metamorphosis of the LC can be 
observed only when the radiative column exists. Hence, 
we may conclude that the existence of the radiative col- 
umn can cause the irregularity of some XBPs. Though it 
is difhcult to explain all the irregularity of XBP LCs by 
only this mixed-beam hypothesis, some of the observed 
complicated LCs may be reconstructed by adopting this 
method. 

5.2. Light curve statistics 

Now, we have a large sample of pulsars listed in Table 1 
- 3. According to Basko, Sunyaev (1976) and the previous 
discussion in this paper, when the X-ray luminosity of 
XBP is larger than L* ^ 5 x 10^^ [erg s~^], the radiative 
column with finite height will appear on the NS surface. 
If such a radiative column with finite height emits both 
the pencil beam and fan beam and we can observe both 
beams, the observed light curve may show an irregular 
behavior. In our sample, however, not all the pulsars with 
Lx > L* are categorized as Irregular group, and a large 
fraction of bright pulsars show Regular LC. In fact we have 
40 bright pulsars with Lx > 5 x 10^^ [erg s~^] and 24 of 
them are Regular members. If all the bright pulsars should 
show the irregular behavior, this fraction is inconsistent. 

Actually, even if the radiative column appears and both 
beams are emitted, it does not mean every pulsar will be 
observed as irregular one. For instance, when the mag- 
netic axis and the line of sight almost coincide, the pulse 
component due to fan beam will be very weak. Hence, the 
observed LC will be dominated by pencil beam only, and 
cannot be irregular. Also when the magnetic axis almost 
coincides with rotation axis, the irregularity will not ap- 
pear since the LC becomes flat. For another example, if 
the beam is coUimated in small angle due to certain rea- 
sons (it corresponds to large fi and in Eq. (4)), we may 
less likely observe both beams in different phase. Also 
in this case, the observed LC will be dominated by only 



(a) (b) 




0.5 1 1.5 2 0.5 1 1.5 2 

phase phase 



(c) (d) 




0.5 1 1.5 2 0.5 1 1.5 

phase phase 



Fig. 6. Model light curves. Pure-fan, pure-pencil, and 
combined beams are shown in the branches entitled fan, 
pencil and total in the figure. The computed parameters 
are the following. Models (a), (b) and (c) are generated 
for 9ni = 7r/8, 7r/4 and Stt/S, respectively. For these three 
cases, the mixing factor is fixed as A = 0.5 and also we fix 
fi = u = I. Through all the figures, solid, dashed, dash- 
dotted, and dotted lines correspond to Bo = 7r/8, 7r/4, Stt/S, 
and n/2, respectively. In model (d), we show the result 
for fi = u = 2 but other values arc the same as model (b). 

one of two beams and behave regularly. Then, the possi- 
bility of observation of Irregular pulsars depends on the 
inclination angles and the opening angles of the beams. 

Suppose that the opening angles of the pencil- and fan- 
beam are ipp and ipf, respectively. When the inclination 
condition 

26i,„ + ¥>f + V3p > I (6) 

is satisfied, there is a region where both the pencil beam 
and the fan beam can be observed from the same position. 
The probability that the observer's location enters in this 
region can be described as 



where 7 = 20^^ + ip[ + ipp — tt/2. p^^ denotes the probability 
that the observer can see both pencil- and fan-beam, and 
in this case the light curve will be categorized as Irregular 
group. Of course, the probability that the light curve is 
observed as a Regular member is Pr ~ 1 — fir. Note that 
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though the beam opening angle depends on the strong 
magnetic field, the gravitational bending, beaming, etc., 
here <p denotes the observable opening angle and contains 
all those effects. 

Next, we consider the probability that n Regular pul- 
sars and m Irregular pulsars are observed, within totally 
(n + m) pulsars. Assuming that the (n + m) sampling 
is random, this probability can be given by the bimodal 
distribution: 



{n + m)\ 



(8) 



n\ ml 

Hence, the probability that the observed number of 
Irregular pulsar is less than M is given by the summa- 
tion 



M 
m— 1 



Pn- 



(9) 



where = n + m is the total sample number. 

In Fig. 7, wc show the probability that we will observe 
M, or less. Irregular pulsars within the total sample N, as 
the function of pir = 2j/tt. When wc consider the bright 
pulsars with Lx > 5 x lO^^, iV = 40 and M = 16. This 
case is shown by thick solid line in the figure. Imposing 
brighter criteria, = 29 and M = 14 for Lx > 1 x 10^"^, 
and A^ = 14 and il/ = 8 for Lx > 5 x 10^'^, respectively. 
Those results are shown by thick dashed and dotted lines 
in the figure. 

Our diagnosis of light curves contains some uncertain- 
ties, and we denote the pulsars whose LC are marginally 
Regular by "Rm" in the tables. If we count those 
marginal pulsars as Irregular members, the data count 
will be changed as M = 26 for Lx > 5 x lO^^, M = 22 for 
Lx > 1 X 10^^, and M = 11 for Lx > 5 x 10^'', respectively. 
The probability functions for those cases arc shown by 
thin lines in the same figure. 

This figure shows that, in order to reconstruct the ob- 
served number count of Irregular pulsars, the probability 
Pir should be small, under the present assumptions. It 
directly means that the summation of the magnetic incli- 
nation angle 9^ and the opening angles of beams iff and 
ifp should be small. In this analysis 9^1 and ip are degen- 
erated. However, in order to explain the large number of 
bright Regular pulsars, anyway magnetic inclination 9^ 
should be small. In order to achieve Pm<M > 0.5, we need 
Pi, < 0.42 {N = 40, M = 16, Lx > 5 x lO^^) for the severest 
case (thick sofid line in Fig. 7), and pi, < 0.82 (A^ = 14, 
M = 11, Lx > 5 x lO^"^, with "Rm"s) for the lax case (thin 
dotted line in the figure). On the other hand, pi, should 
be larger than 0.02 - 0.04, so highly coUimated beams are 
also unacceptable. 

About the magnetic inclination angles, Bulik et al. 
(2003) have analysed the ratio of single peak LCs over 
all the observed LCs, and they conclude that average 6'm 
should be small in order to explain the observed ratio of 
single peak LC. It is interesting and to be emphasised that, 
although our analysis is completely independent from the 
light curve analysis given by Bulik et al. (2003) which is 



0.01 



0.001 




Fig. 7. Probability of observing M Irregular pulsars within 
N samples, as the function of pi^. Thick lines (A, B and C) 
correspond to the XBPs of Irregular group, while thin lines (a, 
b and c) show the Irregular and marginal LC pulsars. Curves 
with marks of (A, a), (B, b) and (C, c) denote results for 
Lx > 5 X lO^"^, Lx > 1 X 10^'^ and Lx > 5 x 10^'', respectively. 

using the pulse peak number, our result also suggests the 
same conclusion of small ^m- 

Another reason why the sample of Regular bright pul- 
sars is over-sized might be the bias of observed energy 
bands. In this analysis, we are assuming that the fan- 
beam and the pencil-beam will be emitted in different en- 
ergy bands. However, the energy bands in our sampling 
may not cover these different energy bands separately. For 
instance, even if assume that, for a certain XBP, we have 
observational data in 2-4keV and 4-lOkeV, it is not always 
that the pencil-beam and the fan-beam could be observed 
separately in those two bands. In our diagnosis, unfortu- 
nately the used energy bands are miscellaneous and are 
not standardized. In actual diagnosis, 9 samples of 24 
bright Regular pulsars are based on the observed LC in 
only 2 energy bands. For the higher reliability, we clearly 
need further samples in broader energy bands from soft 
X-ray to gamma-ray. For bright pulsars, the quality of 
the LC becomes better and small features in the LC can 
be distinguished, hence those bright pulsars may tend to 
be judged as Irregular type. Also from this observational 
bias, this unexpectedly large fraction of Regular compo- 
nent looks a kind of new mystery of XBP LCs. 

In fact, the radiation process from the radiative column 
itself has been still unclear. For the effective analysis, 
we need to know the exact energy band where the X-ray 
photons are emitted for both the pencil- and fan-beam. 
In order to know the reliable predominance of the emis- 
sions due to pencil-beam and fan-beam, we need to solve 
radiative-magneto- hydro equations. Especially, from the 
point of view of beam opening angles, the radiation emit- 
ted from the radiative region in the strong magnetic field 
becomes important. 

6. Summary and conclusion 

In this study, we analyze the observed light curves of 
X-ray binary pulsars. We mainly focus on those XBP sys- 
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terns whose LCs are observed in multiple energy bands. 
We found that some of them show different shapes of 
LC in different energy bands, and wc named such LCs as 
"Irregular" type. Looking for the origin of such an irreg- 
ularity, we found that the X-ray luminosities of Irregular 
systems are higher than ^ 5 x 10'^^[erg s~^]. This crite- 
rion coincides with the critical limit for the existence of 
the accretion column with finite height, given by Basko, 
Sunyaev (1976). The coincidence of the X-ray luminosity 
limit of Irregular XBP with the apparent condition of the 
accretion column is fairly consistent because of following 
reason: When the accretion column exists, the X-ray ra- 
diation will be emitted not only as the pencil-beam, but 
also as the fan-beam. Hence, if we observe those differ- 
ent beams (and also combined beams) in different energy 
bands, the observed X-ray LC will be different in each 
energy band. 

Additionally, we analysed the fraction of Irregular XBP 
in all the observed system. It is apparent that not all 
the XBP with large Lx show Irregular LC. This occurs 
because for the observation of Irregular LCs, the loca- 
tion of the observer has to satisfy the condition given in 
Eq. (6). The probability to take such a location depends 
on the inclination angle of magnetic axis regarding to the 
rotational axis, and the opening angles of beams. To ex- 
plain the observational fraction of Irregular systems in our 
sample, the small inclination angle of the magnetic axis 
is strongly required. This result is consistent with the in- 
dependent analysis given by Bulik et al. (2003). In our 
sampling, however, the completeness of the observed en- 
ergy dependence is not so good. In near future, we would 
like to collect larger number of samples with broad en- 
ergy band, and analyze the data in detail. At the same 
time, an improvement of the theoretical modelling of X- 
ray emission mechanisms in XBP is strongly required for 
better understanding of these systems. 

The author would like to thank Michela Mapelli for her 
helpful discussions, careful reading, and comments. He 
thanks Monica Colpi for her comments. He thanks an 
anonymous referee for careful reading and valuable sug- 
gestions. 
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Table 2. List of X-ray pulsars in SMC 
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Table 3. List of X-ray pulsars in LMC 



No. Name LC Ps Po 

LI RXJ0502.9-6626 — 4.1 — 

L2 LMC X-4 IR 13.7 24.5 

L3 RXJ0529.8-6556 Rm 69.5 — 

L4 SAXJ0544. 1-710 Rm 96.08 — 

L5 EXO054011-65512 — 272 — 

L6 A0535-668 — 0.068 16.7 



4 X 1037 

1 X 1037 

2 X 1036 
2 X 1036 

3.9 X 1035 

1 X 1039 



Bp Comments 

— Schmidtke et al. 1995 

— Woo et al. 1996, Paul et al. 2002, 
Haberl et al. 2003 

— Haberl et al. 1997, 2003 

— Cusumano et al. 1998 

— Haberl et al. 2003 

— Skinner et al. 1982 



